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Semi-global weak stabilization of bilinear 
Schrodinger equations 

Karine Beauchard* Vahagn Nersesyan^'* 



Abstract 



We consider a linear Schrodinger equation, on a bounded domain, with 
biUnear control, representing a quantum particle in an electric field (the 
Qi^ ' control). Recently, Nersesyan proposed explicit feedback laws and proved 

■^C , the existence of a sequence of times (in)neN for which the values of the 

solution of the closed loop system converge weakly in H^ to the ground 
state. Here, we prove the convergence of the whole solution, as i — >■ +00. 
C^ , The proof relies on control Lyapunov functions and an adaptation of the 

LaSalle invariance principle to PDEs. 



Resume: Stabilisation faible semi-globale d'equations de Schrodinger 

Cn ' bilineaires. 

^ , Nous considerons une equation de Schrodinger lincaire, sur un domaine 

borne, avec un controle bilineaire, modelisant une particule quantique 
dans un champ electrique (la commande). Recemment, Nersesyan a pro- 
pose des lois de retroaction explicites et demontre I'existence d'une suite de 
temps {tn)nev> auxquels les valeurs de la solution du systeme boucle con- 

in ' vergent faiblement dans H^ vers I'etat fondamental. Ici, nous demontrons 

la convergence de toute la solution, quand t — >■ +00. La preuve repose sur 
des fonctions de Lyapunov et une adaptation du principe d'invariance de 
LaSalle aux EDP. 
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Version francaise abregee: On considere le systeme 

iz = -Az + V{x)z + u{t)Q{x)z, x e D, (0.1) 

z\aD = 0, (0.2) 

z{0,x) = zo{x), (0.3) 

ou i? C M™ est un domaine borne a bord lisse, V,Q £ C°°(Z3,]R) sont des 
fonctions donnees, u est le controle et z est I'etat. II modelise une particule 
quantique dans un potentiel V et un champ electrique u. 

Notons {ek,v)kGN* les vecteurs propres de I'operateur (—A + V), (—A + 
y)sk,v — ^k.v^k,v, Piyz := z — {z,eiy)eiy la projection orthogonale de 
L'^{D,C) sur Wcct{eky',k > 2} et S' la sphere unite de L'^{D,C)- 

Lcs lois do retroaction explicites suivantes sont introduites dans |15| 

u{z) := -5Im [(a(-A + V)Piy{Qz), (-A + V)Piyz) - {Qz, eiy){eiy, z)] , 

(0.4) 
on S,a > 0. Elles permettent de considerer le systeme boucle 

iz^ -Az + V{x)z + u{z)Q{x)z, x e D. (0.5) 

Rappelons la condition suivante, introduite egalement dans |15] . 

Condition 0.1. Les fonctions V,Q G C°°{D,'R.) verifient 

(i) {Qeiy, ej,v) 7^ pour tout j > 2, 

(ii) Ai V — Aj y 7^ Ap y — Ag y pour tout j, p,q>l tels que {l,j} ^ {p, q} et 



Dans [15], Nersesyan dcmontre que, sous la Condition lO.il il existe une suite 
de temps (tn)nGN auxquels la solution du systeme boucle converge faiblement 
dans H^ vers I'etat fondamental: z{t„) -^ eiy dans H^ quand n -^ +c>o. Dans 
cet article, nous demontrons que toute la solution converge: z{t) -^ ei.y dans 
H^ quand t — >• +00. 



Theorem 0.2. On suppose la Condition W. 1\ veriUee. SoitlAt la resolvante du 
systeme boucle I10.5\) . i0.2\) . II existe un ensemble fini ou denombrable J C K!|. 
tel que, pour tout zq G SDHq DH"^ n'appartenant pas C := {ceiy : c G C, |c| = 
1}, il existe a* := q:*(||zo||2) > tel que, pour tout a 6 (0,a*) — J, Ut{zo) -^ C, 
dans H^ , quand t — > 00. 

La preuve du Theoreme 10.21 repose sur le principe d'invariance de LaSalle 
et se fait en deux etapes. Dans un premier temps, on vcrifie que renscmble 
invariant coincide localement avec C. Dans un dcuxicme temps, on dcmontre la 
convergence. Pour cela, on montre que lcs sculcs valcurs d'adhcrence possibles, 
pour la topologie faible H^, de la solution du systeme boucle, sont dans C. 
Considerant une valcur d'adherence faible H^, Ut^{zo) -^ Zoo, on dcmontre 
qu'elle appartient a C, en montrant qu'elle engendre une solution invariante, 
u[l/lt{zoo)] = 0. Pour cela, on demontre que u[lJt„+tizo)] -^ quand n — >■ +00 
pour presque tout t e [0, +00) et on justifie le passage a la limite [n — )• +cxd] 
dans le feedback. 



1 Introduction 

1.1 The system 

We consider the system (|0.ip - (|0.3[) where D C M™ is a bounded domain with 
smooth boundary, V,Q & C°°{D,M.) are given functions, u is the control, and z 
is the state. It represents a quantum particle in a potential V^ in an electric field 
u. The following proposition establishes the well-posedncss of system (|0.ip - (|0.3p 
(see [5] for a proof). 

Proposition 1.1. For any zq G Hq n H^ (resp. zq G L^) and for any u £ 

L,^Qj,([0, c«),M) problem h0.1\) - W^) has a unique solution z £ C([0,oo),7Jg D 
H^) (resp. z E C{[0,oo), L'^)). Furthermore, the resolving operator Ut{-,u) : 
L" — )■ L" taking zq to z{t) satisfies the relation 

\\Utizo,u)\\ = \\zol Vt>0. (1.1) 

In all this article, ||.|j (resp. ||.||s) denotes the usual norm on L^{D, C) (resp. 
H'{D,C), for every s € W). S is the L2(D,C)-sphere and 



if, 9)-= / f{x)gix)dx,yf,g€L\D,C). 



D 

1.2 Bibliography 

We refer to [2j |4l [71 [M] for exact or approximate controllability results for the 
system (|0.ip - ()0.3p . with open loop controls. This article is concerned with closed 
loop controls: we search explicit feedback laws, that asymptotically stabilize the 
ground state. 

In [13], the same question is addressed for ODE models. The control design 
relies on control Lyapunov functions, and the convergence proof relies on the 
LaSalle invariance principle. This reference deals with the situation where the 
linearized system around the ground state is controllable. The degenerate case 
is studied in [3]. 

The goal of this article is to adapt the result of [T5j to PDE models. In- 
deed, the LaSalle invariance principle is a powerful tool to prove the asymptotic 
stability of an equilibrium for a finite dimensional dynamic system. However, 
using it for infinite dimensional systems is more difficult (because closed and 
bounded subsets are not compact). 

A first possible adaptation consists in proving approximate convergence re- 
sults, as for example in [51112). A second possible adaptation consists in proving 
a weak convergence, as, for example, in [1] and in this article. A third possible 
adaptation consists in proving a strong convergence, as for example in [8]. In 
this case, one needs an additional compactness property for the trajectories of 
the closed loop system. Another strategy consists in designing strict Lyapunov 
functions, as for example in P]. 



1.3 Stabilization strategy 

Let us recall the stabilization strategy proposed in [T5]. We introduce the Lya- 
punov function 

V{z) := a||(-A + V)Pi,vzf + l-\{z, ei,y)p, zeSr\Hlf\ H^, 

where a > 0, {ek,v)k&N' are the eigenvectors of the operator —A + V, (—A + 
V)eky = ^k.v^k,v and Piyz := z — {z, eiy)ei,v is the orthogonal projection 
in L^ onto the closure of Spa.n{eky,k > 2}. Notice that V{z) > for all 
z G S n Hi n H"^ and V{z) = if and only if z = ceiy, \c\ = 1. For any 
z € SfMlln H'^, we have 

V{z) > a\\{-A + V)P,yz\\' > |||A(Pi,v2)|p - Ci > ^\\Azf - C^, 



where Ci and C2 arc positive constants. Thus 

C(l + V(z))>||z||2 



fl.2) 



for some constant C > 0. Following the ideas of [Sj [15], we wish to choose a 
feedback law u{-) such that 

^V(z(t)) < 

for the solution z(t) of ([UH - IJU^)) . Let us assume that Az(t) e H^ Ci H'^ for all 
t > 0. Using dnH]), we get 



-V(z(t)) = 2a Re [((-A + V)Piyz, (-A + V)Piyz) 



2 Re 



{z,eiy){eiy,z) 



2a Re 
-2 Re 



((-A + V)Piy{iAz - iVz - iuQz), (-A + V)Piyz) 
{iAz — iVz — iuQz, eiy){eiy, z) 



Integrating by parts and using the facts that V is real valued, Piy commutes 
with —A + V and 



(-A + V)Piyz\dD = z\dD = eiylao = 0, 



we obtain 



2a Re 



(-^(-A + VfPiyz, (-A + V)Piyz) 



-2 Re 



(lAz — iVz, eiy){eiy, z) 



2a Re 

+ 2aRe 



-lV(-A + V^)Pi,yZ, V(-A + V)Piyz) 
{-lV{-A + V^)Pi,yZ, (-A + V)Piyz) 



2Ai^yRc 



{iz,eiy){eiy,z) 



0. 



Thus 



-V{z{t)) = 2ulm a((-A + V)Piy{Qz), (-A + V^)Pi.vz) - {Qz, eiy){eiy,z) 
at I 



Let us take u{z) defined by (|0.4p where S > 0. Then 

±V{zit)) = -^u\z{t)), (1.3) 

thus t I— >■ V{z(t)) is not increasmg and one may expect that z(t) -^ C := 
{ceiy,c e C, \c\ = 1}, in some sense, when t — > +c«. We consider the closed 
loop system (|0.5p . The following proposition ensures the well posedness of this 
system and the validity of the computations performed above. 

Proposition 1.2. For any zq e H^ n H^ problem W.5\) . W.2\} . Ii0.3\) has a 
unique solution z G C{[0,oo), Hq D H^)- Moreover if Azq G Hq n H^ , then, 
Az€C{[0,oo),H^nH^). 



The local well-posedness and the regularity of the solution of (jO.Sp . (j0.2|) . 
(|0.3p is standard (see [B]). From the construction of the feedback law u it 
follows that a finite-time blow-up in H^ n H^ is impossible. Hence the solution 
is global in time. 

1.4 Main result 

Let us introduce the following condition on the functions V and Q. 
Condition 1.3. The functions V,Q e C°°(D,R) are such that: 
(i) {Qeiy, Cjy) ^ for all j > 2, 
(a) Ai y — Aj y 7^ Ap y — Aq y for all j,p, q > I such that {1, j} ^ {p, q} and 

See the papers [161 El E] for the proof of genericity of this condition. The 
below theorem is the main result of this article 



Theorem 1.4. LetUt be the resolving operator of the closed loop system W. 5\) . 
W.2\) . Under Condition ] 1.3[ there is a finite or countable set J C M!J_ such that 
for any a ^ J and zq £ S H Hq D H^ with < V{zq) < 1 we have 

Ut{zo) ^ C inH'^ast^oo, (1.4) 

where C := {cei_y : c G C, |c| = 1}. 

Remark 1.5. This Theorem proves the semi-global stabilization of the ground 
state. Indeed, for every zq € S Ci Hq D H^ such that zq ^ C, one may chose 
a = a(||2;o||2) > small enough so that the condition < V{zq) < 1 is fulfilled. 

2 Convergence proof 

The first step of the proof consists in checking that the LaSalle invariance set 
locally coincides with C. 



Proposition 2.1. We assume Condition \1.3[ There exists a finite or countable 
set J C M+ such that, for every a ^ J, for every zq G SDHq n H^ with 
{zo, eiy) 7^ and u{Ut{zo)) = for all t > 0, then zq G C. 

This proposition is proved in jl5) . The second step of the proof consists in 
proving the convergence. First, we need the foUowing preliminary resuh. 

Proposition 2.2. Let lAt be the resolving operator of the closed loop system 
W.5\) . W.S]) . Let Zn S Hq D H^ be such that Zn -^ z^o in H^ and Zn -^ Zoo in 
Hq. For every T > 0, there exists N C (0,r) with zero Lebesgue measure such 
that 

1. Ut{zn) -^ Ut{zoo) tn H^ and Ut{zn) ^ Ut{zoo) m H^, Vt e (0,T) - N, 

2. u[Ut{zn)\ ^ u[Ut{z^)], Vt e (0,T) - N. 

Proof of Proposition [2T2l 

First step: Let us show that, if z^ £ Hq D H'^ , z„ -^ z^ in H^ and Zn — > 2:00 in 
Hq, then u{zn) — > u(zoo)- Then, the second conclusion of Proposition \^7I\ will 
be a consequence of the first one. 

Notice that (|0.4p and the fact that Q is real valued imply that 

u{z) = -dlTJi {aQ{-A + V)z,{-A + V)z) +u{z)^u{z), 

where 

u{z) = -Sim [(a(-A + V)PiyiQz), (-A + V){-{z, e,y)eiy)) 

+ {a{-A + V){-{Qz, ei,v)ei,v), (-A + V)z) 

+ {a{-\/Q-Vz-zAQ),{-A + V)z) ^^■^> 

-{Qz,eiy){eiy,z) . 

Thus, passing to the limit in the previous equality, we get u(z„) — > u{zoo). 

Second step: Let us prove the first conclusion of Proposition \2.2\ Let Zn G Hq H 
H^ be such that z„ — ^ z^o in H^ and z„ ^ z^o in H^. For T > define the 
Banach space W := {z ^ C{[{),T],Hl n iJ^) such that z e L^{[0,T], L"^)} en- 
dowed with the norm ||z||w := \\z\\c{[o,T],Hinm) + U\\m[o,T],L^)- The sequence 
of functions (i G [0,T] >-^ iYf(z„))„gpj is bounded in W, and the embedding 
W -^ L^{{0,T),H^) is compact, by Theorem 5.1 in [TO]. Let Y G L^{{0,T),H^) 
and (p an extraction such that 

U.iz^in)) -^Y{.) in L\iO,T),H^). 

Thanks to the Lebesgue reciprocal theorem, one may assume that 

l^tiz^in)) ^ Y{t) in H^yt G (0,r) - N, (2.2) 

where N C (0,r) has zero Lebesgue measure (otherwise take another extrac- 
tion). 



For every t* £ (0,T) — A^, the sequence (Wf«(z^(„)))„gN is bounded in H^ 
and its only possible weak H^ limit is Y{t*) because of (|2.2p . Thus the whole 
sequence converges: V(t'{z^{n)) ~^ Y{t*) in H^. Therefore, we have 

Utiz^^n)) ^ Yit) in i/2 and i^t(z^(„)) ^ y(t) in i/i,Vt G (0,r) - iV. (2.3) 

We deduce from the first step that 

u[Utiz^(^n))] ^ u[Y{t)],yt e (0,T) - N. 

Let A := -A + V^. We fix i* G (0, T) - iV. For every n, we have 

Jo 

Passing to the limit [n —^ +00] in Hq in this equality, using the dominated 
convergence theorem and the continuity of e~ z with respect to z in TJq norm, 
we get 

Jo 
Thus, y(i) = Ut{zao)yt e (Oj'^) ^ ^ (uniqueness of the solution of the closed 
loop system). 

This proves that the sequence {t G [0,T] 1-^ Ut{zn))neH has a unique adher- 
ence value in L^((0,T), TJg). Therefore the whole sequence converges (i.e. one 
may take <^ = Id). We deduce from (|2.3p that 

Ut{zn) -^ Ut{zoc) in H^ and Wt(z„) ^ Ut{zoc) in i/o\Vi ^ (0-^) - ^• 

n 

Remark 2.3. The key point of this proof is that the feedback law u(z) is well 
defined for z strictly less regular that H^ (see (|2.ip : formally, z e H^^^ is 
sufficient) . 

Proof of Theorem II. 4t Let J be as in Proposition 12.11 and a ^ J. Let 

Zo <E S n Hq n H'^ with < V(zo) < 1. Let us prove that the weak H^ w-hmit 
set of {Ut{zo); t ^ 0} is contained in C. 

Let Zoo S Hq n H^ and t„ — >■ +00 be such that Ut,Xzo) ^ ^00 in H^ ■ Let us 
show that Zoo G C. One may assume that Ut,^{zQ) — !■ Zqo in i/g- 

There exists an extraction ip and a subset A^i C (0, +00) with zero Lebesgue 
measure such that 

u[Wt^(„,+t] ^ 0,Vt e (0,+^) - A^i. 

Indeed, the sequence of functions {t £ (0, +00) i-J> u[Ut„+t{zn)])7ieN tends to zero 
in L^(0, +00) because 1 1-^ u[Ut{zQ)] belongs to -L^(0, +00) (see (|1.3p '). 

Let T e [0, +00). Thanks to Proposition H^ there exists N C (0,r) with 
zero Lebesgue measure such that 

M[Wt^,„)+t(zo)] ^ u[Ut{zoo)],yt e (o,r) -iv. 



The uniqueness of the hmit ensures that u[Vlt{zoo)] = 0, Vi G (0, T) — [N Li Ni] . 
Finally, the function t i— >■ u[Ut{zoo)] is continuous and vanishes on (0, T) — 
[NUNi], thus it vanishes on [0, T]. This holds for every T > 0, thus u[llt{zao)] = 
0,Vt G [0, +oo). As V{zoc) < V{zq) < 1, we have (zoo,ei,y) ^ 0. Thanks to 
Proposition 12. 11 we get z^o G C D 

3 Conclusion, open problems, perspectives 

We have proposed explicit feedback laws, that asymptotically stabilize the ground 
state, for the system (|0.ip - (|0.3p . To design the feedback laws, we have used con- 
trol Lyapunov functions. The convergence holds semi-globally in H^ and for the 
weak iJ^-topolopy. The proof relies on a adaptation of the LaSalle invariance 
principle to PDEs. 

Generalizations with different regularities are possible: with Lyapunov func- 
tions inspired by the H'^ distance to the target, one may prove weak H" stabi- 
lization. Generalization to other bilinear equations (for instance wave equations) 
are possible. 

Our proof uses compact injections between Sobolev spaces on a bounded 
domain. Thus, the stabilization question when such compact injections cannot 
be used is still an open problem. 

Another open problem concerns the simultaneous stabilization of N identical 
Schrodinger equations, around N different eigenstates, with only one closed loop 
control. Indeed, if we design feedback laws in the same way as in this article, 
then, the LaSalle invariance set does not coincide with the target. Thus, new 
ideas need to be introduced to tackle this problem. 

The same question for non linear Schrodinger equations is also an open 
problem. 
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